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Abstract 

It is shown that the two-level correlation function R(s, s') in the invariant 
random matrix ensembles (RME) with soft confinement exhibits a "ghost 
peak" at s ~ —s'. This lifts the sum rule prohibition for the level number 
variance to have a Poisson-like term var(n) = rjn that is typical of RME 
with broken U(N) symmetry. Thus we conclude that the U(N) invariance is 
broken spontaneously in the RME with soft confinement, rj playing the role 
of an order-parameter. 
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The statistical description of complex systems by ensembles of random matrices turned 
out to be a powerful general approach that was successively applied to a great variety of 
systems in different fields from nuclear physics [1] to mesoscopics [2] and quantum chaos [3]. 

The classical random matrix theory (RMT) by Wigner, Dyson and Mehta [1] describes 
the statistics of eigenvalues for a Gaussian ensemble of random Hermitian iV x iV-matrices 
H with the probability distribution P(H) oc exp[— TrH 2 ]. By definition, the statistical 
properties of this ensemble are invariant under unitary transformations U(N) of matrices 
H and thus there is no basis preference in the RMT. This means that the classical RMT 
can be applied only to quantum systems where all (normalized) linear combinations of 
eigenfunctions have similar properties. 

For disordered electronic systems, it implies that all eigenstates must be extended. In 
other words, the classical RMT is applicable only for describing the energy level statistics in 
the metal phase [4,5] that exists in the dimensionality d > 2 at a relatively weak disorder. 

With disorder increasing the system goes through the Anderson transition to an insulat- 
ing phase in which all eigenstates are localized. The level statistics in this phase obviously 
cannot be described by the C/(7V)-invariant RME, since one can construct an extended state 
by a linear combination of localized states randomly positioned throughout the sample. 
Thus the proper probability distribution -P(H) must contain a basis preference in order to 
exclude unitary transformations which would lead to formation of such extended states. 

The ensemble of random banded matrices (RBME) [6,7] is an example of such a non- 
invariant RMT. It describes properties of systems belonging to so-called quasi-ld universality 
class which includes quasi-ld disordered electronic systems with localization [6] and certain 
quantum chaotic systems [7]. The corresponding eigenvalue statistics are Poissonic in the 
N — > oo limit (at a fixed bandwidth b) and reduce to the Wigner-Dyson form at b — > oo. 
Thus changing the parameter b/N, one can describe the crossover from Wigner-Dyson to 
Poisson level statistics which occurs in quasi-ld disordered systems with increasing the ratio 
L/£ of the sample length L and the localization radius £. 

While the localization in quasi-ld systems seems to be well described in terms of the 
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RBME, the problem of the random matrix description of the critical region near the Ander- 
son transition and the Anderson insulator phase for d > 2 remains open. The recent works 
[8-10] where the existence of the universal critical level statistics has been demonstrated, 
resulted in an intensive search for the proper RME description. 

In this connection, two different generalizations of the classical RMT have been recently 
proposed [11,12]. The generalized RME studied in Ref.fll] was obtained from the Gaussian 
invariant ensemble by introducing a symmetry-breaking term: 

P(H) oc e- Trii2 e -WTK[A,H][A,H]t)_ (1) 

The /i-dependent term breaks the U(N) invariance and tends to align H with a symmetry 
breaking unitary matrix A thus setting the basis preference. It turned out [11] that even 
after averaging over A the resulting ensemble leads to the eigenvalue statistics that deviate 
from the Wigner-Dyson form. The difference between the Wigner-Dyson statistics that 
correspond to h — and the level statistics for any non zero h turns out to be dramatic 
in the thermodynamic (TD) limit N — > oo. Namely, for h ^ the variance var(n) of the 
number of levels in an energy window that contains n levels on the average, grows linearly 
with n at n > 1: 

var(n) = ((5n) 2 } = rj(h)n ~ hn, r/(0) = 0. (2) 

For the classical RMT [1], var(n) oc Inn that is negligible as compared to Eq.(|2]) for any 
nonzero < t](h) < 1 in the limit n — > oo. The Poisson-like behavior described by Eq.(|2]) is 
valid also for RBME in the TD limit. 

In contrast to Eq.(Jl]), the probability distribution suggested in Ref.[12] is explicitly U (AO- 
invariant: 

P(H)oce- T ^ H ). (3) 

The only singularity in this model is that the (even in E) "confining potential" V(E) grows 
extremely slowly: 
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V(E) = jln 2 \E\, \E\^oo. (4) 

However, the Poisson-like behavior, Eq.(^j), turns out to be valid for this model too, provided 
that the energy window does not contain the origin E = 0. 

It should be stressed [13] that for steeper confining potentials, V(E) = \E\ a , no deviation 
from the Wigner-Dyson statistics was observed in the bulk of the spectrum, so that i] a = 
for all a > 0. What happens with the invariant RME at the transition from a power-law to 
logarithmic confinement, looks like the spontaneous breaking of the {7(iV)-syminetry. 

In this Letter we present both analytical and numerical arguments showing this novel 
phenomenon to exist. We show the parameter 77 to play the role of an order-parameter and 
clarify its connection with the breaking of the normalization sum rule in the TD limit. 

The reason why rj = for a wide class of invariant RME is connected with the normal- 
ization sum rule: 

[ +0 °Yf\s,s>)ds' = l. (5) 

J — oo 

Here the cluster function Y 2 (s, s') = S(s — s') — R(s, s') is related to the two-level correlation 
function R(s, s'): 

(6) 

1 ; (p(E s )) (p{E s ,)) 1 ; 

where p{E) = Tr{S(E — H)} is the level density, (...) denotes ensemble averaging with the 
probability distribution -P(H), and the new variable s(E) is chosen so that the mean level 
density in this variable, (p(s)) = (p(E s ))dE s /ds = 1 for —N/2 < s < N/2. 

For iV finite, Eq.(|5|) is an exact property of any RME that follows simply from the 
normalization condition J p(E)dE = N. However, the sum rule may be violated [14] in 
the limit iV —>■ oo, since after the integration in Eq.([Sp both terms in Eq.(|^) result in the 
divergent constants iV to be subtracted from each other. Below we will assume this limit to 
be taken. 

The general relationship between the cluster function Y^s,*') and the coefficient rj = 
lim n _ +00 {<i[var(n)]/(in} in the level number variance, Eq([|), for an energy window centered 



at a point E = E(sq), reads: 



r\ = lim 



a+ 



rv°° i 



L*2 



a+,s 



') + 



a_, s 



')] ds' 



(?) 



where a± = sq ± n/2. In the case where the cluster function in the TD limit Y£°(s, s') 
Y£°(s — s') is translationally invariant, the relationship reduces to: 



/+oo 
Y£°{s)ds. 
-oo 



(8) 



Now comparing Eq.(||) and Eq.® we can conclude that the parameter r\ is equal to the 
deficiency of the sum rule in the N —>■ oo limit [9]. 

Whether the sum rule is broken or not depends on the behavior of the effective level 
interaction at large distances. Quite generally, the joint probability distribution P[{xj}] of 
eigenvalues Xi of a matrix H can be represented in the form similar to the Gibbs distribu- 
tion P[{xj}] = exp[— flTC] for the one- dimensional plasma of classical particles described by 
the potential energy functional = ^y(ij) + W[{xj}]. For the unitary ensembles 

considered in this Letter, the effective temperature (3 = 2. In case of the invariant ensembles 
given by Eq.(^), the many-body interaction VTfjxj}] reduces to the pair- wise logarithmic 
effective level interaction [1]: 
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(9) 



i>3 



For RME similar to that given by Eq. ([!]), the corresponding effective level interaction that 
appears after averaging [11] over the symmetry-breaking matrix A, contains all the many- 
body terms and can be rewritten in the form of a determinant Detij, i,j = 1, ...N: 



W[{xi}\ = -llnDet^e-^ 2 ^-^]. 



(10) 



In both limiting I — > and I the determinant in Eq.fllOD can be 

approximated by a pair-wise interaction of the form: 



W[{ Xi \] 



i£m(l-e 

i>3 



- 2h2N2 ^-^ 2 )+ const. 



(11) 



The main difference between Eq.([|) and Eg. (pTTj) is that the symmetry-breaking field 
h ytz results in the cut-off of the effective level interaction at sufficiently large distances 
[15]. This difference is crucial for the fulfillment of the sum rule, Eq.(|5]), in the TD limit. 

It can be shown that the normalization sum rule persists also in the TD limit if the 
effective level interaction is long-range (non-integrable). We will present here a simple 
proof of this statement based on the mean-field approach [16] that is valid exactly for the 
long-range interactions in the N — > oo limit. Suppose that the effective level interaction 
f(xi — Xj) is pair-wise and long-range. Then using the relationship 5{p(E))/6V(E r ) = 
-P[(p{E)p{E')) - (p(E)) (p{E'))} and the Dyson mean-field equation [1] for (p(E)) 

/+oo 
{p{E"))f{E" - E') + V(E') = const, (12) 
-oo 

one can derive the integral equation for the two-level correlation function [16]: 

/+oo 
ds"R°°(s, s") f{E s „ - E s ,) = p-His - s'). (13) 
-oo 

Now integrate this equation over all s changing the order of integration in the left-hand side 
and denote the integral R°°(s, s")ds = I: 

/+oo 
f{E s „ - E a ) ds" = 1. (14) 
-oo 

From the definition of the function E s just after Eq.(|6|) and a physically obvious property 
d{p{E))/d\E\ < that holds for an even confinement potential V(E), it follows that E s 
increases linearly or faster with s. Therefore for any long-range interaction with / f(E)dE 
divergent at infinity, the integral in Eq.([14]) is also divergent. Thus we arrive at the rela- 
tionship: 

/+oo 
Y™(s,s')ds' = 0, (15) 
-oo 

which proves the validity of Eq.(|5]) in the limit N — ► oo. On the contrary, if the effective 
interaction f(E-E') is cut at large distances, the quantity I in Eqs.(|l4]),(|T5[) must be finite 
and thus the sum rule is violated. 
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We see that the symmetry-breaking term in Eq.(|l]) leads to the cut-off of the effective 
level interaction at large distances and thus to the break-down of the normalization sum rule 
in the TD limit. This results in a quasi-Poisson level number variance, Eq.(|2|), governed by 
the nonzero parameter r], Eq.([8]). 

The situation resembles the appearance of the long-range order in spin systems in the 
external magnetic field h which breaks down the rotational invariance. In this case, the 
spin-spin correlator (S(r)S(r')) — > m 2 is constant at large distances, the magnetization |m| 
depending linearly on |h|. Since the spin-spin correlator is invariant under global rotations 
of the spins, the order parameter, that is \m\, is unchanged after averaging over the direction 
of the symmetry-breaking field h. 

In our problem, the quantity analogous to (S(r)S(r')), which is invariant under global 
(independent on E, E') unitary transformations, is the two-level correlation function 
(Tr{5(E — H)}Tr{5(E' — H)}). The symmetry-breaking field hA plays the role of the 
magnetic field, and the parameter rj is similar to the magnetization. Like \m\, rj remains 
nonzero after averaging over A. 

Now we return to consider the TD limit of the invariant ensemble with soft confinement, 
defined by Eqs.(|3]),(f|). Since the effective level interaction, Eq.([5|), is long-range, the sum 
rule, Eq.(||D, must hold for any confining potential V(E). However, for soft confinement 
this does not necessarily lead to rj = 0. The reason is the dramatic break-down of the 
translational invariance of the cluster function Y£°(s,s') = Y£°(s — s') + Y a (s,s') which 
turns out to have an anomalous part Y a (s,s f ) that exhibits a sharp peak of width A near 
s = -s'. 

With the translational invariance broken in such a way, the level number variance var(n) 
becomes dependent on the position of the energy window Sq, and can oscillate as a function 
of n. The parameter rj is given by the general Eq.(^), where n = 1,2,.... It is easy to 
show that the integral in Eq.(^) reduces approximately to that in the sum rule, Eq.(||), only 
provided that n/2 — \sq\ ^> A, that is the origin E = s = is far within the energy window. 
This means that the parameter 77 vanishes for a symmetric window sq = and it is nonzero 
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if | s 1 > n/2. 

Indeed, the Monte Carlo simulations on the classical Id plasma described by Eqs.(|]),(|9]) 
show a dramatic difference in the level number variance in these two cases. A remarkable 
property of the model is that for a symmetric window, the level number variance var(ra) is 
constant for all integers n ^> 1. Thus the "level rigidity" is even higher than that for the 
classical RMT where var(n) oc Inn. 

Though the anomalous part Y a (s,s') can be extracted from the exact solution [12] for 
the particular confining potential obeying Eq.(f|), its existence and importance was never 
mentioned before. Below we present a simplified derivation of the cluster function Y£°(s, s') 
for small values of the parameter q = e _7r2j4 . 

We start with the representation [1] of Y£°(s, s') in terms of the orthonormal "wave func- 
tions" <fii{E) = pi(E)exp(—V(E)), where Pi{E) are orthogonal polynomials corresponding 
to the weight function e~ 2V ^: 

2 1 ' J K(E s ,E s )K(E s ,,E s ,y 
Here the kernel K(E, E') is given by [1]: 

K{E ' E) = ^ W^E ' (17) 

where C is a constant and <p (e) are N — > oo limits of wave functions ipn(E) of the odd and 
even order, respectively. 

For q 1, the wave functions <p ( e \(E) have a "semi-classical" form that is a generaliza- 
tion of the iV — ■> oo limit of the oscillator wave functions for a nonlinear s(E): 

<p (E) = Csm[7rs(E)\, <p e (E) = C cos[tts(E)}. (18) 

The arguments in sin and cos are chosen so that the density of zeros po = ds/dE coincides 
with the mean level density (p(E)). 

The function s(E) can be easily found from the known solution of the Dyson mean field 
equation (H^). Solving Eq.flT!^) with f(E — E') = — In \E — E'\ for the confining potential, 



Eq.(|J), and taking the limit iV — > oo, we have: 
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(p(E)) = ds/dE = jfa, E s = \e 2 ^ A sign(s) J (19) 

where A is a constant of integration. 

Now using Eqs.(|16|)-(|T9|) and the relationship [1] (p{E)) = K(E,E) one can obtain an 
explicit form of the cluster function ^(s, s') for q = e~ n2A <C 1, where A is a coefficient in 
Eq.(|J). For ss' > we arrive at the translationally invariant expression [12]: 

Y°°(s -J)-— Sin2[?r(g ~ s/)] ss' > (20) 

However, there are strong correlations for ss' < too: 

^'>^£SfeSr ss ' <0 - (21) 

This is just the anomalous part of the cluster function discussed above. Its remarkable 
property is a sharp peak at s ~ —s' with a height that does not decrease when \s — s'\ ~ 
2|s| — > oo. 

The existence of such a "ghost" peak in the two- level correlation function is confirmed 
also by the Monte Carlo simulations on the effective plasma model with a soft confinement 
[see Fig.l]. Such anomalous, infinitely long-range level correlations make it possible for 
the invariant random matrix ensemble defined by Eqs.([3|), (§]) to exhibit the Poisson-like 
level number variance var(n) = rjn which is typical for RME with broken U(N) symmetry. 
Moreover, after the substitution ir 2 A = 2/h ^> 1, the form of the "normal", translationally 
invariant part of Y£°(s,s') given by Eq.(^) is identical to that found in Ref. [11] for the 
matrix ensemble, Eq. ([!]), with symmetry breaking. 

This is a crucial point in the chain of arguments in favour of the statement that the 
U(N) symmetry is spontaneously broken in the random matrix ensemble, Eqs.(|3]),(|4]), with 
soft confinement. 

The corresponding order parameter which arises either due to the external symmetry 
breaking field h or spontaneously for a sufficiently soft confinement, can be defined in the 
N — ► oo as follows: 
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V 



lim 1 



I, 



s>>0 



F 2 °°(s,s') ds' 



(22) 



It coincides with the coefficient in the quasi- Poisson term in the level number variance, 



In case of the spontaneously broken U (N) symmetry, a complementary order parameter 
can be introduced: 



which involves integration only over the anomalous part of the cluster function. Because of 
the sum rule, Eq.(|5]), the sum rj + fj = must be zero. Exploiting the spin analogy discussed 
above, one can say that the spontaneous symmetry breaking is of the AFM type, with r\ 
and fj playing a role of the sublattice magnetizations. 

We conclude that the anomalous statistics of eigenvalues in RME with soft confinement 
is a result of the spontaneous breaking of the U(N) symmetry which exhibits itself in the 
quasi- Poisson term in the level number variance and in the appearance of the "ghost" peak 
at s — —s' in the two-level correlation function R(s, s'). 
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Eq.(0) provided that the energy window does not contain the origin. 



n 




(23) 
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FIGURES 

FIG. 1. MC results showing the existence of the "ghost" correlation peak. The simulations are 
performed with N = 100 "particles" in the confining potential V(E) = ^ ln 2 (l + |I?|) for (3 = 2, and 
A=0.5. The reference particle was mobile around s = 24. The solid line in the inset corresponds 
to Eq.(21) with oscillations being averaged out. 
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